Wilson versus Clover fermions: A case for improvement 
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We present evidence for improvement with tadpole improved clover fermions based on an analysis of the chiral 
behavior of Bk and the quark condensate. Also presented are a comparison of the mass splittings in the baryon 
octet and decuplet, a calculation of ca using standard 2-point correlation functions, and the problem of zero 
modes of the Dirac operator. 



1. B K 

With Wilson fermions, straightforward calcu- 
lations of Bk using the 1-loop improved AS = 2 
operator fail due to the large mixing with the 
wrong chirality operators [1]. Since this mixing 
is an artifact of lattice discretization, one hopes 
that it can be significantly reduced by improving 
the action. By comparing results obtained using 
the Wilson and the tadpole improved clover ac- 
tion (csw = 1-4785) on the same quenched gauge 
lattices (170 lattices of size 32 3 x 64 at (3 = 6.0) 
we show that this is indeed the case. 

Fig. 1 shows the Wilson and clover data as a 
function of a 2 M K . For each data set, Bk is writ- 
ten as the sum of two parts the contribution 

of the diagonal (the 1-loop tadpole improved LL) 
operator, and the mixing term which is propor- 
tional to a s . The general form, ignoring chiral 
logarithms and terms proportional to (m s — md) 2 , 
for pi = p f = M K is [1] 
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Figure 1. Data for Bk for Wilson and clover 
fermions. The diagonal (AS = 2) and mixing 
operators contribution arc shown separately. Ver- 
tical lines are at M K = 495 MeV. 
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The coefficients a, @, Si are pure artifacts, there- 
fore their value can be used to quantify improve- 
ment. Of these a is the most serious as it causes 
Bk to diverge in the chiral limit. 

The divergence, in the limit Mk — ► 0, of 
the diagonal term due to a non-zero a is evi- 
dent in Fig. 1 for Wilson fermions. This arti- 
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fact is only partially cancelled by the 1-loop mix- 
ing operator. The situation is considerably im- 
proved with clover fermions. The corresponding 
values at M K = 495 MeV are B K (Wilson) = 
-0.29(1) whereas B K (clover) = 0.50(1). This im- 
provement arises because the two dominant arti- 
facts a and (3 are significantly reduced; a(W) — 
0.0460(14) versus a(C) = -0.0125(6), and 
f3(W) = 0.048(46) versus j3(C) = -0.006(37). 

As explained in [1], the contributions propor- 
tional to a, (3, Si can be removed completely by 
studying the momentum dependence of the ma- 
trix elements. Short of calculating the mixing co- 
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efficients non-perturbatively, the way to remove 
the artifacts in 7, 62, S3 is to extrapolate to a = 0. 
We have done the calculation at f3 = 6.0 only, 
where our final results are Bk{ NDR, 2 GeV) = 
0.68(7) and 0.58(7) for Wilson and clover formu- 
lations respectively. The benchmark value, in- 
cluding a — > extrapolation, is Bk (Staggered) = 
0.628(41), as obtained by the JLQCD collabora- 
tion [2]. 

2. QUARK CONDENSATE 

The chiral condensate (iptp) is not simply re- 
lated to the trace of the Wilson quark propaga- 
tor (Sf(0, 0)). The breaking of chiral symme- 
try by the r term introduces contact terms that 
need to be subtracted non-perturbatively from 
(S^ (0,0)) [3]. This has not proven practical. In- 
stead, the methods of choice are to either evaluate 
the right hand side of the continuum Ward Iden- 
tity 

(^/>} WI = (0|S F (0,0)|0) 

= lim m q [ d i x(0\P(x)P(0)\0), (1) 
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or cast the Gell-Mann, Oakes, Renner relation 
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in terms of lattice correlation functions [4] . These 
estimates have errors of both 0(a) and 0(ma), 
and at fixed a are therefore expected to agree only 
in the chiral limit. A comparison of the efficacy 
of the two methods is shown in Fig. 2. 

We find that a reliable extrapolation to the chi- 
ral limit can be made using a linear fit, and the 
two methods give consistent results for both Wil- 
son and clover fermions. Also, the 0(ma) correc- 
tions are significantly smaller for clover fermion. 

3. BARYON MASS-SPLITTINGS 

In Ref. [5] we presented a detailed analysis 
of mass-splittings in the baryon octet and de- 
cuplet with Wilson fermions. We had found a 
large non-linear dependence on quark mass for 
the £ — A, E — N, and 3 — N splittings. Extrap- 
olation of the data to the physical masses includ- 
ing these non-linearities gave estimates consistent 
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Figure 2. Data for (i/jip) versus the quark mass for 
Wilson (W) and clover (C) fermions. The vertical 
lines show m u ,d and m s for the two actions. The 
l/a are 2.33(4) and 1.99(5) GeV for Wilson and 
Clover fermions. 

with observed values. On the other hand we had 
found a surprisingly good linear fit to the decuplet 
masses, and the splittings were underestimated 
by 25 — 30%. The data with clover fermions show 
the same qualitative features. As an illustration, 
we show a comparison of the £ — A splitting in 
Fig. 3. Details of the analysis will be published 
elsewhere [6]. 

4. DETERMINATION OF c A 

The improvement coefficient for the axial cur- 
rent, ca, is calculated using the the axial WI [7]. 
If the clover coefficient csw is tuned to its non- 
perturbative value 1.769 at (3 = 6.0 [7], the sum 
(mi + TO2) of quark masses defined by 

£ ^ [A M + ac A d^P]W (x, t)jW (0)) 

E^ (12) (<W (21) (o)> 

should be independent of t and the initial pseu- 
doscalar state created by J, up to corrections of 
0(a 2 ). We vary the composition of the initial 
state by using J — P or A4 and by using "wall" 
or "Wuppertal" smearing functions in the calcu- 
lation of the quark propagators. The results in 
Fig. 4 show a large dependence on the initial state 
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Figure 3. Behavior of the E — A splitting versus 
the sum of the three quark masses. The leading 
mass dependence m s — (m u + md)/2 has been 
divided out. 

for Wilson fermions and almost none already for 
c sw = 1-4785! We estimate c A = -0.026(2) from 
this clover data, whereas the ALPHA collabora- 
tion report c A = -0.083(5) at c SW = 1-769 [7]. 
We are repeating the calculation at csw — 1.769 
to understand this difference. 

5. ZERO MODES 

The explicit breaking of chiral symmetry in 
Wilson-like fermions gives rise to the problem 
of "exceptional configurations" in the quenched 
theory. The cause is that the Wilson r term 
breaks the anti-hcrmitian property of the mass- 
less Dirac operator. As a result, zero modes of 
the Dirac operator extend into the physical region 
k < k c . Thus, on a given configuration, as the 
quark mass is lowered and approaches the first of 
the unphysical modes, one encounters exception- 
ally large fluctuations in the correlation functions. 
Such configurations dominate the ensemble aver- 
age and as discussed in [8] there is no basis for ex- 
cluding them. Tuning csw reduces the 0{a) chi- 
ral symmetry breaking artifacts as shown above, 
however, it does not reduce this problem [8]. We 
find, by comparing fluctuations in 2-point and 3- 
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Figure 4. Wilson and clover fermion results for 
ca as a function of (mi + 771.2), and for three dif- 
ferent sources J. The clover data is seen to be 
independent of J and shows a linear behavior. 

point correlation functions between Wilson and 
Clover fermions, that the problem, in fact, gets 
worse. A deeper understanding of the persistence 
of the zero mode problem even though the chiral 
behavior is improved is missing. 
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